o 
o 

(N 
o 

Q 

in 



Oh! 



Multiplicity fluctuations in hadron-hadron and 
nucleus-nucleus collisions and percolation of 

strings 

Pedro Brogueira*and Jorge Dias de Deus*,^ 
February 2, 2008 



> 
O 

CN ■ Abstract 
<N , 

CN I We argue that recent NA49 results on multiparticle distributions 

^ ■ and fluctuations, as a function of the number of participant nucleons, 

■ suggest that percolation plays an important role in particle production 

^ I at high densities. 



Recentely, the NA49 collaboration has presented results, from the CERN/SPS 



Qj [ at 158 A GeV, on multiplicity fluctuations or, to be more precise, on V{n)/ < 



n >, 



X ■ V{n)/ < n <-'>-<->\ (1) 

H : ^ " <n> ^ ^ 

as a function of the number Npart. of participant nucleons, from pp to PbPb 
collisions [1]. 

These data are very interesting for several reasons: 

1) They show evidence for universal behaviour: the experimental points 
in the plot V{n)/ < n> versus Npart. fall into a unique curve (see Fig.l). 

2) The physics in the small Npart. limit {pp, Npart. 2) and in the large 
Npart. limit {PbPb, Npart. — > '^'ApbPb) seems to be quite the same, as in both 
cases the quantity (1) approaches 1. The fluctuations are larger in the inter- 
mediate Npart. region (see Fig.l). 
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3) The (negative) particle distribution, in the low density and in the high 
density limits, is in fact a Poisson distribution (see Fig. 2), the distribution 
being wider than Poisson in the intermediate Npart. region. 

In the framework of the string model with percolation [2], these results 
are quite natural. On one hand, percolation is a universal geometrical phe- 
nomenon, the properties depending essentially on the space dimension (di- 
mension 2, impact parameter plane, in our case), and being controlled by the 
transverse density variable 77, 

n^ijps, (2) 

where r is the transverse radius of the string (r ~ 0.2/m), R the radius of 
the interaction area, and Ns the average number of strings. The quantity 
{R/rY is nothing but the interaction area in units of the string transverse 
area. As R and Ns are functions of the number Npart. of participants, Npart., 
similarly to f], becomes, at a given energy,a universal variable. 

On the other hand, in percolation [3], what matters is the fluctuation in 
the size of the clusters of strings: one starts, at low density (small Npart.), 
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from a situation where strings are isolated, at intermediate density one finds 
clusters of different sizes, and one ends up, at high density, above the perco- 
lation threshold, with a single large cluster. In both, low and high, density 
limits, fluctuations in cluster size vanish (see Fig.3). In the simplest string 
model the particle distribution is Poisson (as observed in e'^e~ and pp at low 
energy) and V/ < n >— > 1 in both, low and high, density hmits (see Figs.l 
and 2). 
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Let us try to be more specific. In hadron-hadron and nucleus-nucleus 
collisions, during the collision strings are produced along the collision axis, 
and these strings may overlap and form clusters of different sizes. In the spirit 
of percolation theory, we shall assume that fluctuations in the number of 
strings per cluster dominate over fluctuations in the number Nc of clusters. 

Generalizing [4], we write for the multiplicity distributions, for a given 
value of r] or N^art.-i 



Pin) = J dNcipiN) J U^^,dNdl^,W{Ni)Pr,,{iii)d (n - E/^^j > (3) 

where (p{Nc) is the probability of having Nc clusters, W{Ni) the probability 
of having a cluster of strings, PN.{iJ,i) the probability of a cluster of 
strings emitting /li particles. We shall further write 



N 



(4) 



where Pi{n) is the single string distribution, assumed to be Poisson. 

Note that in (3) we treat the cluster as independent, and in (4) we treat 
the strings themselves as independent (this is to be corrected later by includ- 
ing string fusion). 

As mentioned above, we neglect fluctuations in the number of clusters. 
This means that we approximate (fi{Nc) by a 5-function, such that (3) be- 
comes 

P{n) = J uf^,dNidl^iW{Ni)Pr,Md (n - E/^^j , (3') 

where Nq is the average number of clusters, for a given density rj, or number 
of participant nucleons Npart. ■ 
Prom (3') and (4) we obtain, 

<n>^ Nc < N >n , (5) 

and 

< > - < n >2= Nc[{< N^>-<N >^)n^+ < N > {n^ - n^)] , (6) 
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where n is the single string particle multiplicity, and, finally, for the (vari- 
ance) / (multiplicity) ratio. 



(8) 



<n> < N > 

where 

V{N) _ <N^> - <N>^ 

<N> ^ < > 
The previous percolation argument for the fiuctuations in the size of the 
clusters (see Fig. 3) requires: 

aud ^ ^ . (9) 



<:N>r]^0 <N> 1] ^ oo 






ri«1 rjssl r|»1 

<N>%1 <N>?^N3 

In percolation there are two sum rules to be satisfied. If N stands for the 
number of strings in a cluster, Nc the number of clusters and Ns the number 
of strings, we have: 

Nc < N >= Ns , (10) 

and 

Na < A >= - e-'^) , (11) 

r 

where < A > is the average area occupied by a cluster (see, for instance, 
[3,5]). In percolation < A > increases with t], approaching the full area of 
interaction in the rj oo limit. 

In the case of percolation we thus write for < A>, [5], 

<A>=fm-ni-e'^)-l] + l, (12) 
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where f{r)) is a percolation function, such that /(ry) — > 0, as — > 0, and 

< A >— > 1, as expected for isolated strings, and — 1, as ?7 — oo, and 

< A >—> (— )^ as expected in percolation. For f{r]) we have chosen 



f{rj) ^ {1 + e-^^-^^^^T' , 



(13) 



with a = 0.85 and 1]^ = 1.15 [5]. 

From (10), (11) and (12) we obtain 

<^>=Y^(/('7)[(7f(l-e-) 

with < >-> 1 as ?7 — > 0, and < N >^ Ns, as t] - 
Regarding the cluster string variance, < A^^ > 
satisfy the general constraint, 

< A^^ > - < >^> , 

and in addition we imposed the constraints, 

< A^^ > - < A^ >^ 

and 



(14) 



oo (see Fig. 3). 
- < N >^, we have to 



(15) 



?7 ^ 







< A^^ > - < A^ >^ 



We thus write for the variance, 
V{N) =<N'^ > - <N 



V 



oo 



. 



1 1 



[1 + bri)e' 



-brj 



<N>' 



(16) 
(17) 

(18) 



where 6 > is an adjustable parameter (fixed at the value b — 1.65). Note 
that (18) satisfies (15), (16) and (17). 

Before making a comparison between our string percolation model and 
NA49 data, there are two questions to be addressed: 

i) F{'r]) factor due to random colour summation 

When strings fuse in a cluster the effective colour charge is not just the 
sum of the colour charges of the individual strings [6] . In pratice, the effective 
number A^ of strings is reduced, [7] , 



N 



'<A> 
<N> 



N — > F{7])N , 



(19) 
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where (see (10) and (11)), 



1 - e-'^ 

Fiv) = , (20) 

such that, instead of (5) and (7), have 

< n F{r])Nc <N>n^ F{7])Nsn , (5') 

and 



<n> ^ ' < N > 

Note that the F{r}) correction is more important for < n >, (5'), then for Eq. 
(7'). If the single particle distribution is Poisson with average multiplicity n, 
the average cluster has also a Poisson distribution with multiplicity F{r)) < 
N > n. 

ii) The relation between rj and Npart. 

In the definition of 77, (2), what appears is not Npart. but rather the average 
number Ns of strings and the radius R of interaction. Making use of simple 
nuclear physics and multiple scattering arguments, one has [8] 

R ~ RpN'f , (21) 

and 

Ns ^ N^N'f , (22) 

where Rp is the nucleon radius (~ Iffn), Ng is the (energy dependent) 
number of strings in pp collisions, at the same energy, and Na is given by 

Na^^. (23) 

We would like to mention that Eqs.(21) and (22) are not rigorous: in (21) 
geometrical factors are not taken into account, in (22) no distinction is made 
between valence strings and sea strings, [8]. 

Prom (2), (21), (22) and (23) we obtain for the relation between rj and 

Npart. , 

^)'^sNT (24) 
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We shall now present our results: 

1- V(n)/ < n > 

In Fig.l we show our curve, Eq.(7') with (14), (18) and (24), in comparison 
with NA49 data. The obtained values for n{n = 0.12) and iVf(A^f = 4.5) are 
consistent with the additional constraint < n >p~ 0.52 (as seen in Fig.2). 
Note that the behaviour of our curve is essentially determined by the be- 
haviour of (< A/"^ > - < A/" >^) / < A/" >. 

In Fig.2 we show fits of the multiplicity distributions for different values 
of Npart., with Negative Binominals. Reasonable fits are obtained with values 
for the NB parameter: k = oo, Poisson, at low and high density, and k — 29, 
for intermediate density. 

2- < n >Na and < n >p 

We can relate < n >Na to < n >p by making use of (5') and (24): 

< n ^ < n >, nT (25) 

we first note that (25) satisfies saturation when A^ oo: 

<n>NA — ' const. (26) 

Relation (25) is valid for high energy. At low energy - ^ ~ 20GeV is the 
energy at SPS - the presence of valence quarks cannot be ignored. We take 
them into account by writing, instead of (25), 



nvNA) 



< " >^-= < ^ V-^i}- ^I^T)] , (250 

r {Up) 

where c is a parameter decreasing with energy, 1 > c > 0, measuring the 
relative contribution of valence quarks to multiplicity, 

c^^, (27) 

< n >p 

and for c was taken the value c = 0.59. In Fig.4 we show (25') in comparison 
with NA49 data. 
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3- Pt fluctuations 

It has been argued [9] that a convinient quantity to measure genuine Pt 
correlations - above random correlations - is 0Pj, 



<n> 



-y z 



(28) 



where, for a given event, Z = Yh=i {^Ti — Pt), z = Pt — Pr, Pt meaning 
averaging over the single particle inclusive distribution. NA49 has presented 
results for [10], in the same conditions of the multiplicity experiment. 
A crude percolation model for Pt fluctuations, developed in [5], gives 



(l>Pr = Pt 



- VjN) 



(29) 



where i — 
in (29) 



Prp Pq 



T — ^ T I /Pt- we compare (28) with data. Note that 



at small and high densities, as seems to be the case in data. 
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In a recent paper [11] an attempt was made to relate 4>Pj.toV{n)/ < n >. 
However a direct proportionality, as proposed in [11], does not seem to satisfy 
the low and high density limits. We would like to mention that the use of 
percolation to study Pt fluctuations was initiated with [12]. 

4- Rapidity interval and centre of mciss energy dependence of 

< n > and k 

The average multiplicity < n >, (5'), 

< n >= F{r])Nc < N > n , 
and the Negative Binomial Parameter k, 

(30) 



or, from (7'), 



< iV >2 



are both proportional to Nq, the average number of clusters. If we assume 
that the clusters of strings are approximately uniformely distributed in the 
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central rapidity region - this is the central rapidity plateau idea - then Nc 
is, for Ay > 0, a linear function of Ay (but approching 1 as Ay 0, 
i.e., < >— 1). The multiplicity is also a linear function of Ay (and 
approaching zero as Ay — > 0, because n is also a linear function of Ay, for 
small Ay). The expected < n > and k dependence on Ay - linear increase 
with Ay for large Ay - agrees with data [13]. 

On the other hand, < n > is an increasing function of y^, because 
Nc < N >= Ns increases with i/i and, at large 1], F{t]) only decreases 

as 1/ \J Ns- Increasing the energy, at fixed number of participants, means 
increasing rj. Regarding k, Nq and < N >^ /(< > — < N >^), these are 
not monoatonic functions of rj. The NB parameter k goes to infinity at small 
T) (energy) and at large r) (energy). One thus expects k, for instance in pp, 
to have a minimum at some energy [14]. 

In conclusion, we find that the recent NA49 results, regarding the mul- 
tiplicity distribution dependence on the number of participant nucleons are 
quite consistent with the impact parameter, percolation description of hadron- 
hadron and nucleus-nucleus coUisions at high energies and high densities. 

We would like to thank Elena Ferrciro, Carlos Pajares and Roberto Ugoc- 
cioni for many discussions, and to thank P. Seybot and M. Rybczynski for 
information on NA49 data. This work has been done under the contract 
POCTI/36291/FIS/2000, Portugal. 
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Figure Captions: 

Fig.l Variance over average multiplicity, for negative particle production, as 
a function of the number of participants. The curve is from (7') with 
(14), (18) and (24). Data are from NA49 [1]. 

Fig.2 Multiplicity Distributions, P(n_), as a function of n_. The curves are 
Poisson (dashed lines) and Negative Binomial (full hue). 

Fig.3 Impact parameter percolation. For small densities [r] <^ 1) and for 
large densities {rj ^ 1) there are no strong fluctuations in the number N 
of strings per cluster. For intermediate densities (?7 ~ 1) A'"-fluctuations 
are large. 

Fig.4 The average multiplicity divided by l/2Npart. as a function of Npart.- 
The curve corresponds to Eq.(25'). Multiphcities were calculated from 
the distributions of Fig.2. 

Fig. 5 (f) (Pt) as a function of Npart.- The curve corresponds to Eq.(29). Data 
are from [10]. 
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